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We formulate the equations of equilibrium of neutron stars taking into account strong, weak, 
electromagnetic, and gravitational interactions within the framework of general relativity. The nu- 
clear interactions are described by the exchange of the cr, cj, and p virtual mesons. The equilibrium 
conditions are given by our recently developed theoretical framework based on the Einstein-Maxwell- 
Thomas- Fermi equations along with the constancy of the general relativistic Fermi energies of par- 
ticles, the "Klein potentials", throughout the configuration. The equations are solved numerically 
in the case of zero temperatures and for selected parameterizations of the nuclear models. The 
solutions lead to a new structure of the star: a positively charged core at supranuclear densities 
surrounded by an electronic distribution of thickness ~ h/{mec) ~ W^h/{m^c) of opposite charge, 
as well as a neutral crust at lower densities. Inside the core there is a Coulomb potential well of 
depth ~ m^c^/e. The constancy of the Klein potentials in the transition from the core to the 
crust, impose the presence of an overcritical electric field ~ {m„/me)^Ec, the critical field being 
Ec = m^(? /{eh). The electron chemical potential and the density decrease, in the boundary in- 
terface, until values < /i™'^'^ and pcrust < Pcore- For each central density, an entire family of 
core-crust interface boundaries and, correspondingly, an entire family of crusts with different mass 
and thickness, exist. The configuration with pcrust = Pdrip ~ 4.3 x 10^^ g/cm'' separates neutron 
stars with and without inner crust. We present here the novel neutron star mass-radius for the 
especial case pcrust = Pdrip and compare and contrast it with the one obtained from the traditional 
Tolman-Oppenheimer-Volkoff treatment. 



I. INTRODUCTION 

It is well known that the classic works of Tolman [l| 
and of Oppenheimer and Volkoff for short TOV, ad- 
dresses the problem of neutron star equilibrium configu- 
rations composed only of neutrons. For the more general 
case when protons and electrons are also considered, in all 
of the scientific literature on neutron stars it is assumed 
that the condition of local charge neutrality applies iden- 
tically to all points of the equilibrium configuration (see 
e.g. [3|). Consequently, the corresponding solutions in 
this more general case of a non-rotating neutron star, 
are systematically obtained also on the base of the TOV 
equations. 

In general, the formulation of the equilibrium of sys- 
tems composed by different particle species must be es- 
tablished within the framework of statistical physics of 
multicomponent systems; see e.g. Q. Thermodynamic 
equilibrium of these systems is warrantied by demand- 
ing the constancy throughout the configuration of the 
generalized chemical potentials, often called "electro- 
chemical" , of each of the components of the system; see 
e.g. 0-01 ■ Such generalized potentials include not only 
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the contribution due to kinetic energy but also the contri- 
bution due to the potential fields, e.g. gravitational and 
electromagnetic potential energies per particle, and in the 
case of rotating stars also the centrifugal potential. For 
such systems in presence of gravitational and Coulomb 
fields, global electric polarization effects at macroscopic 
scales occur. The balance of the gravitational and electric 
forces acting on ions and electrons in ideal electron-ion 
plasma leading to the occurrence of gravito-polarization 
was pointed out in the classic work of S. Rosseland [sj. 

If one turns to consider the gravito-polarization effects 
in neutron stars, the corresponding theoretical treatment 
acquires remarkable conceptual and theoretical com- 
plexity, since it must be necessarily formulated consis- 
tently within the Einstein-Maxwell system of equations. 
O. Klein, in [5|, first introduced the constancy of the gen- 
eral relativistic chemical potential of particles, hereafter 
"Klein potentials" , in the study of the thermodynamic 
equilibrium of a self-gravitating one-component fiuid of 
neutral particles throughout the configuration within the 
framework of general relativity. The extension of the 
Klein's work to the case of neutral multicomponent de- 
generate fluids can be found in I6j and to the case of 
multi-component degenerate fiuid of charged particles in 
0. 

Using the concept of Klein potentials, we have re- 
cently proved the impossibility of imposing the condition 
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of local charge neutrality in the simplest case of a self- 
gravitating system of degenerate neutrons, protons and 
electrons in /3-equilibrium it has been shown that 
the consistent treatment of the above system implies the 
solution of the general relativistic Thomas-Fermi equa- 
tions, coupled with the Einstein-Maxwell ones, being the 
TOV equations thus superseded. 

We have recently formulated the theory of a system of 
neutrons, protons and electrons fulfilling stron g, e lectro- 
magnetic, weak and gravitational interactions [10|. The 
role of the Klein first integrals has been again evidenced 
and their theoretical formulation in the Einstein-Maxwell 
background and in the most general case of finite temper- 
ature has been there presented, generalizing the previous 
results for the "non- interacting" case f£|. The strong in- 
teractions, modeled by a relativistic nuclear theory, are 
there described by the introduction of the cr, uj and p 
virtual mesons |lll - [l4| (see Subsec. FlI Al for details). 

In this article we construct for the first time the equi- 
librium configurations of non-rotating neutron stars fol- 
lowing the new approach, 0, [l^]. The full set of the 
Einstein-Maxwell-Thomas-Fermi equations is solved nu- 
merically for zero temperatures and for selected param- 
eterizations of the nuclear model. We use units with 
7i = c = 1 throughout the article. 



II. THE CONSTITUTIVE RELATIVISTIC 
EQUATIONS 

A. Core Equations 

It has been clearly recognized that, since neutron stars 
cores may reach density of order ^ 10^^-10^'' g/cm"^, 
much larger than the nuclear density pnuc ~ 2.7 x 10^^ 
g/cm'^, approaches for the nuclear interaction between 
nucleons based on phenomenological potentials and non- 
relativistic many-body theories become inapplicable (see 
[3 ^ self-consistent relativistic and well-tested 

model for the nuclear interactions has been formulated 
Within this model the nucleons interact with 
cr, Lo and p mesons through Yukawa-like couplings and 
assuming flat spacetime the equation of state of nuclear 
matter has been determined. However, it has been clearly 
stated in [l^] that, when we turn into a neutron star 
configuration at nuclear and supranuclear, the global de- 
scription of the Einstein-Maxwell-Thomas-Fermi equa- 
tions is mandatory. Associated to this system of equa- 
tions there is a sophisticated eigenvalue problem, espe- 
cially the one for the general relativistic Thomas-Fermi 
equation is necessary in order to fulfill the global charge 
neutrality of the system and to consistently describe the 
confinement of the ultrarelativistic electrons. 

The strong interactions between nucleons are described 
by the exchange of three virtual mesons: a is an isoscalar 
meson field providing the attractive long-range part of 
the nuclear force; w is a massive vector field that models 
the repulsive short range and; p is the massive isovector 



field that takes account surface as well as isospin effects 
of nuclei (see also [H, HBl)- 

The total Lagrangian density of the system is given by 



C — Cg + Cf + C„ + + £p + + £int, 
where the Lagrangian densities for the free-fields are 

C = 
' 16^G' 

r — 5_ p pi^" 

^ ~ 16^ ' 



(1) 

(2) 
(3) 
(4) 
(5) 



i^p = -\n^.n^'' + \mlp,p'^, (6) 

where il^^ = d^uj^, - d^uj^, TZ^^ = 9pp^ - d^p^, = 
df_iA^ — dyA^ are the field strength tensors for the w'', p 
and A^^ fields respectively, stands for covariant deriva- 
tive and R is the Ricci scalar. We adopt the Lorenz 
gauge for the fields ^1^, w^, and p^. The self-interaction 
scalar field potential U{a) is a quartic-order polynom for 
a renormalizable theory (see e.g. 

The Lagrangian density for the three fermion species 

is 



i=e,N 



(7) 



where ipN is the nucleon isospin doublet, ipe is the elec- 
tronic singlet, rrii states for the mass of each particle- 
specie and Dfj^ = 9^ + F^, being F^ the Dirac spin con- 
nections. 

The interacting part of the Lagrangian density is, in 
the minimal coupling assumption, given by 



where the conserved currents are 



1,N 



N- 



(8) 



(9) 
(10) 

(11) 
(12) 



The coupling constants of the cr, w and p-fields are g^, 
and gp, and e is the fundamental electric charge. The 
Dirac matrices 7'^ and the isospin Pauli matrices satisfy 
the Dirac algebra in curved spacetime (see e.g. [l^ for 
details). 

We first introduce the non-rotating spherically sym- 
metric spacetime metric 

ds" = e'^^^Ue - e^('^)dr2 - r^dS^ - sin^ ed^\ (13) 
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where the z/(r) and A(r) are only functions of the radial 
coordinate r. 

For very large number of fermions, we adopt the mean- 
field approximation in which fermion-fleld operators are 
replaced by their expectation values (see 19] for details). 
Within this approximation, the full system of general rel- 
ativistic equations can be written in the form 
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constant, 



(14) 
(15) 

(16) 
(17) 
(18) 



(19) 
(20) 
(21) 
(22) 



where we have introduced the notation cjo — i^, Po — P, 
and Aq = V for the temporal components of the meson- 
fields. Here pi — dS/dni = V {PfY + and rii = 
{Pf)^/{3TT'^) are the free-chemical potential and number 
density of the i-specie with Fermi momentum Pf. The 
particle effective mass is rfiN — m^r + and rrie — 
and the effective potentials Vp^n are given by 



9pP 
9pP- 



(23) 
(24) 



The constancy of the generalized Fermi energies E^ ^ 
Ep and E^ , the Klein potentials, derives from the ther- 
modynamic equilibrium conditions given by the statis- 
tical physics of multicomponent systems, applied to a 
system of degenerate neutrons, protons, and electrons 
within the framework of general relativity (see [10] for 
details). These constants are linked by the /3-equilibrium 
between the matter constituents 



Ei 



(25) 



The electron density is, via Eq 

g-3./2 



Tie 



37r2 



20|) . given by 
[l>2 + 2meF-m2(e'^-l)]3/2, (26) 



where V = eV + Ef . Substituting Eq.([26l) into Eq. ^ 
one obtains the general relativistic extension of the rel- 
ativistic Thomas-Fermi equation recently introduced for 
the study of compressed atoms [l^, [2l| . This system of 
equations has to be solved with the boundary condition 
of global neutrality; see 0, [lo| and below for details. 

The scalar density n^, within the mean- field approxi- 
mation, is given by the following expectation value 



= {''pNtpN) = 



(27r)3 ^ 

i—n.p 



(27) 



where ei{p) — y/p^ + rhi is the single particle energy. 

In the static case, only the temporal components of 
the covariant currents survive, i.e. {il)(x)^^il){x)) = 0. 
Thus, by taking the expectation values of Eqs. (l9))- (IT2l) . 
we obtain the non- vanishing components of the currents 

Jo'' = nchUa = [rip - ne)uo, (28) 
Jq ^ HbUQ ^ {rin + np)uo, (29) 
Jq ^a'^o = {np - nn)uo, (30) 

where nt — rip + rin is the baryon number density and 
''^o — \/goo = 6^/2 is the covariant temporal component 
of the four- velocity of the fluid, which satisfies u'^u^ = 1. 

The metric function A is related to the mass M (r) and 
the electric field E{r) = -c^t'^+^^/^l/' through 

2GM{r) 



-A(r) ^ ^_ 



1 - 



r 

2GM{r) 



Gr^E^{r) 
GQ-'ir) 



(31) 



being Q{r) the conserved charge, related to the electric 
field by Q{r) = r^E{r). 

The energy-momentum tensor of free-fields and free- 
fermions T^'' of the system is 



" 7 



where 



rpflU 

7 



P 



1 

47r 



pv 



(32) 

(33) 
(34) 

(35) 

(36) 
(37) 



where the energy-density £ and the pressure V are given 
by 



V 



E 



(38) 
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being Ei and Vi the single fermion fluid contributions 
2 /"^^ 

= 7^ / e«(p)4Vrfp, (39) 

1 2 f-^^ d2 
= / -7-T47rp2dp. (40) 

3 (27r)3 7o e,(p) 

It is worth to recall that the equation of state (I38l) - p0l) 
satisfies the thermodynamic law 

£+V= '^'Z^- (41) 

The parameters of the nuclear model, namely the cou- 
pling constants Qs, and gp, and the meson masses 
ma-, rriaj and rup are usually fixed by fitting experimen- 
tal properties of nuclei, e.g. saturation density, binding 
energy per nucleon (or experimental masses), symmetry 
energy, surface energy, and nuclear incompressibility. In 
TableUwe present selected fits of the nuclear parameters. 
In particular, we show the following parameter sets: NL3 
[H, NL-SH 1231, TMl [H, and TM2 [H. 





NL3 


NL-SH 


TMl 


TM2 


(MeV) 


508.194 


526.059 


511.198 


526.443 


(MeV) 


782.501 


783.000 


783.000 


783.000 


rUp (MeV) 


763.000 


763.000 


770.000 


770.000 


9s 


10.2170 


10.4440 


10.0289 


11.4694 




12.8680 


12.9450 


12.6139 


14.6377 


9p 


4.4740 


4.3830 


4.6322 


4.6783 


92 (fm-i) 


-10.4310 


-6.9099 


-7.2325 


-4.4440 


93 


-28.8850 


-15.8337 


0.6183 


4.6076 


C3 


0.0000 


0.0000 


71.3075 


84.5318 



TABLE I: Selected parameter sets of the a-uj-p model. 

The constants 32 a-nd 53 are the third and fourth order 
constants of the self-scalar interaction as given by the 
scalar self-interaction potential 

U{cy) = im>2 + lg2^3 ^ 1^^^4 (42) 

The non-zero constant C3 that appears in the TMl and 
TM2 models corresponds to the self-coupling constant 
of the non-linear vector self-coupling \cz{uj ^oj'^Y . We 
have not include such a self-coupling vector interaction 
in the general formulation presented above. However, we 
show also here the results of the integration when such 
a self-interaction is taken into account and we refer to 
[2^ [is'l for details about the motivations of including 
that contribution. 

The numerical integration of the core equations can be 
started given a central density and the regularity con- 
ditions at the origin; see below Sec. IIIII for details. At 
nuclear density the phase-transition to the "solid" crust 
takes place. Thus, the radius of the core i?coie is given by 



£{r = i?coro)/c^ = Pnuc- Thcsc equations must be solved 
with the boundary conditions given by the fulfillment of 
the condition of global charge neutrality and the continu- 
ity of the Klein potentials of particles between the core 
and the crust. 



B. Core-crust transition layer equations 

In the core-crust interface, the mean-field approxima- 
tion for the meson-fields is not valid any longer and 
thus a full numerical integration of the meson-field equa- 
tions of motion, taking into account all gradient terms, 
must be performed. We expect the core-crust transition 
boundary-layer to be a region with characteristic length 
scale of the order of the electron Compton wavelength 
'-^ Ae = h/{mec) ~ 100 fm corresponding to the electron 
screening scale. Then, in the core-crust transition layer, 
the system of equations P^ - (P^ reduces to 
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2 

f -V 
r 




(43) 
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= e^™" [daU{a)+gsns] , 


(44) 
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= -e^"^"" [gujJo - mluj] , 


(45) 


P ^ 


2 , 

- -P 
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= -e^'^"'^ [gpj(l - mlp] , 


(46) 


Ai« = 


= Aip H 


— eV = constant , 

-I- eV + guiUj + gpp — constant , 

-Aie + 2<?ppe-''— 


(47) 
(48) 
(49) 



due to the fact that the metric functions are essentially 
constant on the core-crust transition layer and thus we 
can take their values at the core-radius e'^'^"" = e''^^'^"™-' 
and e-^"^"" = e^''^''°"=K 

The system of equations of the transition layer has a 
stiff nature due to the existence of two different scale 
lengths. The first one is associated with the nuclear 
interactions ^ = ?i/(m^c) ~ 1.5 fm and the sec- 
ond one is due to the aforementioned screening length 
^ — h/irriec) ~ 100 fm. Thus, the numerical inte- 
gration of Eqs. P5t - P5|) has been performed subdivid- 
ing the core-crust transition layer in the following three 
regions: (I) a mean-field-like region where all the fields 
vary slowly with length scale ^ Ae, (II) a strongly inter- 
acting region of scale ~ A^ where the surface tension due 
to nuclear interactions dominate producing a sudden de- 
crease of the proton and the neutron densities and, (III) 
a Thomas-Fermi-like region of scale ^ Ae where only a 
layer of opposite charge made of electrons is present pro- 
ducing the total screening of the positively charged core. 
The results of the numerical integration of the equilib- 
rium equations are shown in Fig. [T][2]for the NL3-model. 

We have integrated numerically Eqs. ((Ti| - (P^ for the 
models listed in Table H] The boundary conditions for 
the numerical integration are fixed through the following 
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procedure. We start assuming a value for the central 
baryon number density 71.^(0) = nn(0) + np{0). From the 
regularity conditions at the origin we have e"^(°) = 1 
and ne(0) — np(0). 

The metric function v at the origin can be chosen ar- 
bitrarily, e.g. 1^(0) = 0, due to the fact that the system of 
equations remain invariant under the shift v — > 1/+ con- 
stant. The right value of i' is obtained once the end of the 
integration of the core has been accomplished and duly 
matched to the crust, by fulfilling the following identity 
at the surface of the neutron star. 



2GM{R) 



(50) 



being M{R) and R the total mass and radius of the star. 
Then, taking into account the above conditions, we solve 
the system (I17p -(f22 |) at the origin for the other unknowns 
(7(0), p(0), n„(0), np(0), ^^(0). 

The initial conditions for the numerical integration of 
the core-crust transition layer equations are determined 
by the final values given by the numerical integration 
of the core equations, i.e. we take the values of all the 
variables at the core-radius i?corc- 

In the region I the effect of the Coulomb interaction 
is clear: on the proton-profile we can see a bump due to 
Coulomb repulsion while the electron-profile decreases as 
expected. Such a Coulomb effect is indirectly felt also by 
the neutrons due to the coupled nature of the system of 
equations. However, the neutron-bump is much smaller 
than the one of protons and it is not appreciable in Fig.[T]- 
[2] due to the plot-scale. In the region II we see clearly 
the effect of the surface tension due to nuclear interac- 
tion which produces a sharp decrease of the neutron and 
proton profiles in a characteristic scale ^ Att. In addi- 
tion, it can be seen a neutron skin effect, analogous to 
the one observed in heavy nuclei, which makes the scale 
of the neutron density falloff slightly larger with respect 
to the proton one, in close analogy to the neutron skin 
effect observed in neutron rich nuclei, see e.g. [1^. The 
region III is characterized by a smooth decreasing of the 
electron density which resembles the behavior of the elec- 
trons surrounding a nucleus in the Thomas-Fermi model. 
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FIG. 1: Upper panel: electric field in the core-crust transi- 
tion layer in units of the critical field Ec. Lower panel: par- 
ticle density profiles in the core-crust boundary interface in 
units of cm~^. Here we use the NL3-model of Table U and 
Act = h/(mac) ~ 0.4 fm denotes the sigma-meson Compton 
wavelength. The density at the edge of the crust in this ex- 
ample is /Ocrust = pdrip = 4.3 X 10^^ g/cHl^. 
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FIG. 2: The same as Fig. [U but setting = in order to see 
the efi'ects of the p-meson with respect to the case gp 7^ 0. 



The matching to the crust must be done at a radius 
^corc + 5R where full charge neutrality of the core is 
reached. Different thicknesses 5R correspond to differ- 
ent electron Fermi energies Ef . The thickness of the 
core-crust transition boundary layer 5R as well as the 
value of the electron density at the edge of the crust, 
^crust _ ne(i?corc + 5R), depend on the nuclear parame- 
ters, especially on the nuclear surface tension. 

The equilibrium conditions given by the constancy of 
the Klein potentials (I^Ujl - ip^ throughout the configura- 
tion, impose in the transition layer the following conti- 
nuity condition 



where = fie{Rcorc), el^™™ = eV{Rco,c), and 

= Me(i?coro + SR), and e^""- ~ e''— . 
In the boundary interface, the electron chemical po- 
tential and the density decrease: /i^''"*'* < /z™"' and 
Pcrust < PcoiG- For each central density, an entire fam- 
ily of core-crust interface boundaries exist each one with 
a specific value of SR: the larger the pcrust, the smaller 
the SR. Correspondingly, an entire family of crusts with 
different mass and thickness, exist. From the continuity 
of the electron Klein potential in the boundary inter- 
face given by Eq. (|5ip. it follows that different values of 
Pcrust ^ correspond to different values of the electron 
Fermi energy > 0. In close analogy to the compressed 
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atoms studied in [l^, the case = corresponds to 
the "free" (uncompressed) configuration, where 6R — ^ oo 
and Pcrust = 0, i.e. a bare core. In this configuration 
the electric field reaches its maximum value. The case 
> is analogous to the one of the compressed atom 
|20| . In Fig. Owe have plotted the electron distribution 
in the core-crust boundary interface for selected densities 
at the edge of the crust pcrust = [pdrip, 10^°, 10^] g/cm^, 
where pdrip ~ 4.3 x 10^^ g/cm^ is the neutron drip den- 
sity. 

The configuration with Pcrust = Pdrip separates neutron 
stars with and without inner crust. In the so-called inner 
crust, the neutrons dripped from the nuclei in the crust 
form a fluid that coexist with the nuclei lattice and the 
degenerate electrons [13]. For definiteness, we present in 
this article the results for configurations Pcrust ^ Pdrip, 
i.e for neutron stars possessing only outer crust. The 
construction of configurations with pcrust > Pdrip needs 
to be studied in more detail and will be the subject of a 
forthcoming work. 
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FIG. 3: Distribution of electrons in the core-crust boundary 
interface for different densities at the edge of the crust, pcrust. 
The larger the pcrust, the smaller the electric field E and the 
smaller the thickness of the interface SR. 



In Figs. [T] and [5J we show the core-crust transition 
layer for the NL3 model of Table U with and without 
the presence of the p-meson respectively. The presence 
of the p-meson is responsible for the nuclear asymmetry 
within this nuclear model. The relevance of the nuclear 
symmetry energy on the structure of nuclei and neutron 
stars is continuously stressed in literature; see e.g. [28l - 
l32j |. The precise value of the nuclear symmetry energy 
plays here a crucial in determining the precise value of the 
p-meson coupling which, in the present case, is essential 
in the determination of the intensity of the electric field 
in the core-crust boundary interface; as can be seen from 
the comparison of Figs. [T] and [5] 



C. Crust equations 

Turning now to the crust, it is clear from our recent 
treatment of white dwarfs [2l| that also this problem can 
be solved by the adoption of Wigner-Seitz cells and from 
the relativistic Feynman-Metropolis- Teller (RFMT) ap- 
proach [2^ it follows that the crust is clearly neutral. 
Thus, the structure equations to be integrated are the 
TOV equations 



dV G{£ +V){M + 4Trr^V) 



dr 
dr 



r2(l 



2GM^ 



(52) 
(53) 



where M = M{r) is the mass enclosed at the radius r. 

The effects of the Coulomb interaction in "solid" -like 
electron-ion systems appears only at the microscopic level 
e.g. Dcbye-Hueckel screening in classical systems ,3^ and 
Thomas- Fermi screening in the degenerate case |34| . In 
order to analyze the effects of the microscopic screening 
on the structure of the configuration we will consider two 
equations of state for the crust: the locally neutral case or 
uniform approximation (see e.g. 351) and, for simplicity, 
instead of using the RFMT EoS 20], we use as second 
EoS the one due to Baym, Pethick and Sutherland (BPS) 
[27| . which is by far the most used equation of state in 
literature for the description of the neutron star crust 
(see e.g. [1]). 

In the uniform approximation, both the degenerate 
electrons and the nucleons distribution are considered 
constant inside each cell of volume V^s- This kind of 
configuration can be obtained only imposing microscop- 
ically the condition of local charge neutrality 



Z 



(54) 



The total pressure of the system is assumed to be en- 
tirely due to the electrons, i.e. 



V = V,= 



3 (27r?i)-' 7o ^Jc^p'^ + mlc^ 



dp, (55) 



while the total energy-density of the system is due to 
the nuclei, i.e. £={A/Z)mj\ini,, where nipf is the nucleon 
mass. 

We turn now to the BPS equation of state. The first 
correction to the uniform model, corresponds to aban- 
don the assumption of the electron-nucleon fluid through 
the so-called "lattice" model which introduces the con- 
cept of Wigner-Seitz cell: each cell of radius i?ws con- 
tains a point-like nucleus of charge +Ze with A nucleons 
surrounded by a uniformly distributed cloud of Z fully- 
degenerate electrons. 

The sequence of the equilibrium nuclides present at 
each density in the BPS equation of state is obtained by 
looking for the nuclear composition that minimizes the 
energy per nucleon for each flxed nuclear composition 
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{Z, A) (see Table H] and ^ for details). The pressure V 
and the energy-density £ of the system are, within this 
model, given by 



E_ _ Wn + Wl ^ Se{nbZ/A) 
Ub A Ub ' 

where the electron energy-density is given by 
2 



(27r)3 



(56) 
(57) 

(58) 



and Wn{A, Z) is the total energy of an isolated nucleus 
given by the semi-empirical formula 

Wn = mnc'^{A- Z)+mpC^Z -bA, (59) 

with b b eing the Myers and Swiatecki binding energy per 
nucleon [36[. The lattice energy per nucleus Wl is given 
by 



Wr = -- 



1.819620Z2e2 



(60) 



where the lattice constant a is related to the nucleon 
density npf by npfU^ = 2. 



III. NEUTRON STAR STRUCTURE 
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FIG. 4: Upper panel: electric field in the core-crust transition 
layer, in units of the critical field Ec- Middle panel: particle 
density profiles in the core-crust boundary interface, in units 
of cm~^. Lower panel: density profile inside a neutron star 
with central density p(0) ~ 5pnuc- We compare and contrast 
the structural differences between the solution obtained from 
the traditional TOV equations (locally neutral case) and the 
globally neutral solution presented here. We use here the NL3 
nuclear parametrization of Tableland Aa = h/{mcrc) ~ 0.4 
fm, denotes the sigma-meson Compton wavelength. In this 
example the density at the edge of the crust is pcrust ~ pdrip = 
4.3 X 10" g/cm^ 



In the traditional TOV treatment the density and the 
pressure are a priori assumed to be continuous as well 
as the local charge neutrality of the system. The dis- 
tinguishing feature of our new solution is that the Klein 
potentials are constant throughout the three regions; the 
core, the crust and the transition interface boundary. An 
overcritical electric field is formed and consequently a 
discontinuity in density is found with a continuous total 
pressure including the surface tension of the boundary. 
In Figs. 13] and [SI we compare and contrast the density 
profiles of configurations obtained from the traditional 
TOV treatment and with the treatment presented here. 

In Figs. [Bl IT^ we show the results of the numerical in- 
tegration of the system of the general relativistic consti- 
tutive equations of the configuration from the center all 
the way up to the surface with the appropriate boundary 
conditions between the involved phases. In particular, we 
have plotted the mass-radius relation as well as the com- 
pactness of the neutron stars obtained with the models 
listed in Table HI 

It is worth to note that the inclusion of the Coulomb 
interaction and in particular the presence of the negative 
lattice energy Wl results in a decreasing of the pressure 
of the cells. Such an effect, as shown in Fig. l9HT2l leads 
to a decreasing of the mass and the thickness of the crust 
with respect to the uniform-approximation case where no 
Coulomb interactions are taken into account. 



Comparing the mass and the thickness of the crust ob- 
tained with these two different EoS, we obtain systemati- 
cally crusts with smaller mass and larger thickness when 
Coulomb interactions are taken into account. This re- 
sults are in line with the recent results in [2l|, where 
the mass-radius relation of white-dwarfs has been calcu- 
lated using an EoS based on the relativistic Feynman- 
Metropolis- Teller model for compressed atoms |2C| . 

In the case of the BPS EoS, the average nuclear com- 
position in the outer crust, namely the average charge 
to mass ratio of nuclei Z/A, is obtained by calculating 
the contribution of each nuclear composition present to 
the mass of the crust. We exemplified the analysis for 
two different cores: Mcore = 2.56M0, i?corc = 12.79 km; 
AfcoiG = 1.35M0, i?coro = 11.76 km. The relative abun- 
dance of each nuclide within the crust of the star can be 
obtained as 

^■^■^TTms I ^^r^Sdr, (61) 

-""crust JAr 

where the integration is carried out in the layer of thick- 
ness Ar where the particular nuclide is present; seeHJand 
Fig. [131 Our results are in agreement with the analysis on 
the neutron star crust composition obtained in (stI . [38| . 
In both cases we obtain as average nuclear composition 
3g''^Br. The corresponding crusts with fixed nuclear com- 
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edge of the crust is pcrust = 10^° g/cm^. 
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the star radius R. In the crust we have used the BPS equation 
of state and the nuclear models are in Table IT] 



IV. OBSERVATIONAL CONSTRAINTS ON THE 
MASS-RADIUS RELATION 



position Is^Br for the two chosen cores are calculated ne- 
glecting Coulomb interactions (i.e. using the first EoS). 
The mass and the thickness of these crusts with fixed 
35^ Br are different with respect to the ones obtained us- 
ing the full BPS EoS, leading to such average nuclear 
composition. For the two selected examples we obtain 
that the mass and the thickness of the crust with aver- 
age Is'^Br are, respectively, 18% larger and 5% smaller 
with respect to the ones obtained with the correspond- 
ing BPS EoS. This result shows how small microscopic 
effects due to the Coulomb interaction in the crust of the 
neutron star leads to quantitative not negligible effects 
on the macroscopic structure of the configuration. 



It has been recently pointed out that the most up- 
to-date stringent constraints to the mass-radius relation 
of neutron stars are provided by the largest mass, the 
largest radius, the highest rotational frequency, and the 
maximum surface gravity, observed for pulsars 39] . 

So far, the highest neutron star mass measured with a 
high level of experimental confidence is the mass of the 
3.15 millisecond pulsar PSR J1614-2230, M = 1.97 ± 
O.O4M0, obtained from the Shapiro time delay and the 
Keplerian orbital parameters of the binary system [iof. 
The fitting of the thermonuclear burst oscillation light 
curves from the accreting millisecond pulsar XTE J1814- 
338 weakly constrain the mass-radius relation impos- 
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TABLE II: pmax is the rtieiximum density at which the nuclide is present;A Tii, A R2 and R.A.1(%), R.A.2(%) are rispectively 
the thickness of the layer where a given nuclide is present and their relative abundances in the outer crust for two different 
cases: Mcore = 2.56M0, i?core = 12.79 km; Mcore = 1.35M0, Rcovc = 11.76 km. 
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FIG. 9: Mass of the crust as a function of the compactness 
for the crust EoS without Coulomb interactions. 



FIG. 10: Crust-thickness as a function of the compactness for 
the crust EoS without Coulomb interactions. 



ing an upper limit to the surface gravity of the neu- 
tron star, GM/{c^R) < 0.24 [4l|. A lower limit of 
the radius of RX J1856-3754, as seen by an observer 
at infinity = R[l - 2GM/{c^R)]-^/'^ > 16.8 km, 
has been obtained from the fit of the optical and X- 
ray spectra of the source [4^; it gives the constraint 
2GM/c^ > R- being i?™ = 16.8 km. As- 

suming a neutron star of M = I.4M0 to fit the Chandra 
data of the low-mass X-ray binary XT, it turns out that 
the radius of the star satisfies R — 14.5^i[ g km, at 90% 
confidence level, corresponding to Roo = [15.64, 18.86] 
km, respectively (see [43| for details) . The maximum ro- 
tation rate of a neutron star taking into account both the 
effects of general relativity and deformations has been 
found to be i/max = lO45(M/M0)i/2(lOkm/i?)3/2 Hz, 



largely independent of the equation of state [42]. The 
fastest observed pulsar is PSR J1748-2246ad with a ro- 
tation frequency of 716 Hz [45|, which results in the con- 
straint M > 0.47(i?/10km)3MQ. In Fig. [Hwe show ah 
these constraints and the mass-radius relation presented 
in this article. 

As discussed by J. E. Triimper in [s^, the above con- 
straints strongly favor stiff equations of state which pro- 
vide high maximum masses for neutron stars. In addi- 
tion, putting all of them together, the radius of a canon- 
ical neutron star of mass M — 1.4M0 is highly con- 
strained to the range i? ^ 12 km disfavoring, at the 
same time, the strange quark hypothesis for these spe- 
cific objects. It is clear from Fig.[T3]that the mass-radius 
relation presented here is consistent with all the observa- 
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V. COMPARISON WITH THE TRADITIONAL 
TOV TREATMENT 



FIG. 12: Crust thickness as a function of the compactness for 
crust with the BPS EoS. 



tion constraints, for all the nuclear parametrizations of 
Table HI We present in Table lllli the radii predicted by 
our mass-radius relation for a canonical neutron star of 
M = \AMq as well as for the millisecond pulsar PSR 
J1614-2230, M = 1.97 ± 0.04Mq. 
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13.82 



TABLE III: Radii (in km) predicted by the nuclear 
parametrizations NL3, NL-Sh, TMl and TM2 of TableHl for a 
canonical neutron star of M = 1.4M0 and for the millisecond 
pulsar PSR J1614-2230, M = 1.97 ± O.O4M0. 



In the traditional TOV treatment local charge neu- 
trality as well as the continuity of the pressure and the 
density in the core-crust transition are assumed. This 
leads to explicit violation of the constancy of the Klein 
potentials throughout the configuration (see e.g. Q). In 
such a case there is a smooth transition from the core to 
the crust without any density discontinuity and therefore 
the density at the edge of the crust is ^ pnuc ^ 2.7 x lO^'' 
g/cm'^. The so-called inner crust in those configurations 
extends in the range of densities /Odrip ^ P ^ Pmic while, 
at densities p < Pdrip, there is the so-called outer crust. 

In Figs. [T5] and [TBI we compare and contrast the mass 
and the thickness of the crust as obtained from the tradi- 
tional TOV treatment with the new configurations pre- 
sented here, discussed previously in Sees. |lT] and |llT] char- 
acterized by Pcrust = Pdrip- 

The markedly differences both in mass and thickness 
of the crusts (see Figs. [TSl and fT6|) obtained from the tra- 
ditional Tolman-Oppenheimer-VolkofF approach and the 
new equilibrium configurations presented here, leads to 
a very different mass-radius relations which we compare 
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FIG. 14: Constraints on the mass-radius relation given by 
J. E. Triimper in [s^l and the theoretical mass-radius relation 
presented in this article in Fig. [G] The solid line is the upper 
limit of the surface gravity of XTE J1814-338, the dotted- 
dashed curve corresponds to the lower limit to the radius of 
RX J1856-3754, the dashed line is the constraint imposed by 
the fastest spinning pulsar PSR J1748-2246ad, and the dotted 
curves are the 90% confidence level contours of constant Roo 
of the neutron star in the low-mass X-ray binary X7. Any 
mass-radius relation should pass through the area delimited 
by the solid, the dashed and the dotted lines and, in addition, 
it must have a maximum mass larger than the mass of PSR 
J1614-2230, M = 1.97 ± O.O4M0. 
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FIG. 15: Mass of the crust given by the traditional locally 
neutral Tolman-Oppenheimer-Volkoff treatment and by the 
new globally neutral equilibrium configurations presented in 
this article. We use here the NL3 nuclear model, see Tabled 



FIG. 16: Thickness of the crust given by the traditional locally 
neutral Tolman-Oppenheimer-Volkoff treatment and by the 
new globally neutral equilibrium configurations presented in 
this article. We use here the NL3 nuclear model, see Tabled 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



Global Neutrality 
Local Neutrality 



» ^ 
1 1 
1 ' 
I ' 
/ ' 
f I 
I I 
/ / 



11 



13 
R (km) 



15 



17 



and contrast in Fig. [T71 



FIG. 17: Mass-Radius relation obtained with the traditional 
locally neutral TOV treatment and with the new globally neu- 
tral equilibrium configurations presented here. We use here 
the NL3 nuclear model, see Table U 



VI. CONCLUDING REMARKS 

We have formulated the equations of equilibrium of 
neutron stars based on our recent works [13] and [1, HO, 
[2l| . The strong, we&k, electromagnetic, and gravitational 
interactions are taken into due account within the frame- 
work of general relativity. In particular, the strong in- 
teractions between nucleons is described by the exchange 
of the (J, uj, and p mesons. The equilibrium conditions 
are given by the set of Einstein-Maxwell-Thomas-Fcrmi 
equations and by the constancy of the general relativistic 
Fermi energies of particles, the Klein potentials, through- 
out the configuration. 

We have solved these equilibrium equations numeri- 
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cally, in the case of zero temperatures, for the nuclear 
parameter sets NL3 NL-SH TMl [H], and TM2 
[23]; see Table H] for details. 

A new structure of the star is found: the positively 
charged core at supranuclear densities is surrounded by 
an electronic distribution of thickness ^ h/(mec) ^ 
10^?i/ (rn^c) of opposite charge and, at lower densities, 
a neutral ordinary crust. 

In the core interior the Coulomb potential well is 
^ m^c^/e and correspondingly the electric field is ^ 
(mj,/mpianck)('TT'7r/'TT'e)^^'c ~ 10~^^Ec. Duc to the equi- 
librium condition given by the constancy of the Klein 
potentials, there is a discontinuity in the density at the 
transition from the core to the crust, and correspondingly 
an overcritical electric field ~ (m^/me)^£'c develops in 
the boundary interface; see Fig. [THU 

The continuity of the Klein potentials at the core-crust 
boundary interface leads to a decreasing of the electron 
chemical potential and density, until values /i"'"^* < /z™"' 
and Pcrust < Pcorc at the edge of the crust, where global 
charge neutrality is achieved. For each central density, 
an entire family of core-crust interface boundaries and, 
correspondingly, an entire family of crusts with different 
mass and thickness, exist. The larger Pcrust, the smaller 
the thickness of the interface, the peak of the electric 
field, and the larger the mass and the thickness of the 
crust. The configuration with pcrust = Pdrip ^ 4.3 x 10^^ 
g/cm'^ separates neutron stars with and without inner 



crust. The neutron stars with pcrust > Pdrip deserve a 
further analysis in order to account for the reduction of 
the nuclear tension at the core-crust transition due to the 
presence of dripped neutrons from the nuclei in the crust. 

All the above new features lead to crusts with masses 
and thickness smaller than the ones obtained from the 
traditional TOV treatment, and we have shown specifi- 
cally neutron stars with pcrust = Pdrip; see Figs. [TSHTOl 
The mass-radius relation obtained in this case have been 
compared and contrasted with the one obtained from the 
locally neutral TOV approach; see Fig. [T71 We have 
shown that our mass-radius relation is in line with ob- 
servations, based on the recent work by J. E. Triimper 
[111; see Fig. [H for details. 

The electromagnetic structure of the neutron star pre- 
sented here is of clear astrophysical relevance. The pro- 
cess of gravitational collapse of a core endowed with elec- 
tromagnetic structure leads to signatures and energet- 
ics markedly different from the ones of a core endowed 
uniquely of gravitational interactions; see e.g. (46l - l49| . 

It is clear that the release of gravitational energy in the 
process of gravitational collapse of the core, following the 
classic work of Gamow and Schoenberg (see [13, l5l| ) , is 
carried away by neutrinos. The additional nuclear and 
electromagnetic energy ~ 10^^ erg of the collapsing core 
introduced in this article are expected to be carried away 
by electron-positron plasma created in the overcritical 
electromagnetic field in the collapsing core. 
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